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I. INTRODUCTION 

The pion transition form factor (TFF) between a photon and 
a pion is extracted from the e + <?~ — » e + e~n° process where 
the n° is produced via the two-photon production mechanism. 
This transition is represented as a function of the photon vir- 
tualities as F 7 fiy t y t (q 2 ,q 2 ). The TFF is then extracted when 
one of the electrons is tagged. This electron emits a highly 
off-shell photon with momentum transfer q\ = —Q 2 and is 
detected while the other, untagged, is scattered at a small an- 
gle and then its momentum transfer q\ is near zero. The pion 
transition form factor is then defined as F K 0y,y*{—Q 2 ,Q) = 

The TFF was measured in the CELLO Q and CLEO Q 
experiments in the momentum transfer ranges 0.7 — 2.2 GeV 2 
and 1.6 - 8.0 GeV 2 , respectively. At 2009, the BABAR Col- 
laboration extended these measurements in the Q 2 range from 
4 to 40 GeV 2 0. And recently 0], the Belle Collaboration 
has measured the form factor in the same BABAR's energy re- 
gion with slightly different results on the high-energy region. 

At low transferred momentum, the TFF can be described by 
the expansion: 

F K « ry {Q 2 ) =a o fl+ a7[ ^+b K ^ + 0(Q 6 )) , (1) 

where the parameter ao can be determined from the axial 
anomaly J5] |6) in the chiral limit of QCD, ao = ^rj- with 
f n the pion decay constant. 

The parameter a K , the slope of the TFF, was measured 
by 0, E, and @, with the results a % = -0.1 1(3) (8), a K = 
0.026 (24) (48), and a K = 0.025 (14) (26) respectively (the first 
error is statistic and the second systematic). The CELLO 
Collaboration estimated a K to be a % — 0.0326(26). wa ,(26) 4ys 
in Ref. [1| using an extrapolation from the region of large 
space-like momentum transfer assuming a vector meson dom- 
inance (VMD) and using at zero transferred momentum the 
current experimental value for the partial decay width Y K o^yy 

(which as we will see later is related to F n 0y*y{Q 2 = 0), 
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the axial anomaly). The KTeV Collaboration also predicted 
a K = 0.040(40) through a model-dependent fit to time-like 
data ITOl . The CELLO prediction however dominates the 
number quoted by the PDG jlT) since the direct measure- 
ments are less precise. 

A VMD fit to all the available data (CELLO, CLEO, 
BABAR and Belle) would yield a K = 0.0275(5) with a 
X 2 jd.o.f . = 2.4 (d.o.f. meaning "degrees of freedom"), 
which means a 1.4 standard deviations from the CELLO re- 
sult. This result suggests that the high-energy data may be 
important for determining low energy properties of the TFF. 

One immediately comes to the question on how to improve 
on the quality of the fits to stabilize the predicted result and 
also on how to assign a systematic error to the fit procedure. 

In Ref. [ 12 1 it was suggested that the VMD is a first step on 
a sequence of particular rational approximations called Pade 
Approximants (PA). In that reference was also suggested that 
using Pade Approximants as a fitting functions to analyze the 
pion vector form factor in the space-like region, one can go 
beyond the VMD in a systematic approximation. 

In the TFF case, this fact is of particular interest since the 
data from the BABAR collaboration cannot be easily accom- 
modated in the VMD picture. With the help of these rational 
approximants, one could reach systematically the intermedi- 
ate and high-energy experimental data producing, at the same 
time, an accurate results for the slope and curvature of the TFF 
at low-energies. 

The Pade technics provides with a simple, model- 
independent and systematic method of fitting data with a 
larger range of convergence than the simple polynomial fit or a 
VMD-like fit (such as the one used by CELLO collaboration 
to extract the a K parameter). Given a function f(z) defined 
in the complex plane, the PA Pj^(z) are ratios of two poly- 
nomials Rn(z) and Qm(z) (with degree N and M resp.) the 
coefficients of those exactly coincides with the coefficients 
of the Taylor expansion of f(z) up to the highest order, i.e., 
f(z)-P^(z) = ^(z N+M+l ). 

The PA method also provides with an estimation of a sys- 
tematic error and could be also used to evaluate the impact of 
the vector excitations in the process considered. The technics 
described here were applied also on the search of resonance 
poles in Ref. |[T3l . 

There are several types of PA but as pointed out in Ref. Ifl2ll 
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the analytic properties of the function to be approximated de- 
termines what PA should be used. The time-like region is 
largely dominated by the p —meson contribution. The natu- 
ral choice seems to be a Pi(Q 2 )- On the other hand, since 
it is well known [ 14 1 that the TFF behaves like 1/Q 2 at very 
large energies one could try to incorporate this information by 
considering a P$ +l (Q 2 ) . 

For explanatory reasons we show here how to constract a 
P\{Q 2 ) approximant [a P$ + i(Q 2 ) is more involved and less 
illustrative] . 

Given a function f(z) defined in the complex plane, a Pade 
Approximant P^ is defined fi"5l . without any loss of general- 
ity, by 

Ptfrzo) = X>(Z-Z0)* + : gfc^ , , (2) 
k=0 1 a L \ Z Zo > 

where the coefficients are the Taylor coefficients of the cor- 
responding f(z) function that is been approximated. 

Eq. Q shows that the pole s p of each P^ is determined by 
the ratio s p = ai/ai + \. 

The TFF seems to be well described by the simple VMD 
ansatz. VMD relies on the accurate knowledge of the light 
meson spectra. When the spectral information is given in 
advance one should also take advantage of that information 
and consider other kinds of rational approximants. These are 
the Pade-type approximants. In the Pade-type approximants 
(PTAs) the poles of the Pade are fixed to certain values (in our 
case, the resonances of the spectrum). 

The simplest PTA sequence incorporates the lowest reso- 
nance, the M p , and it is called Tj . The famous VMD ansatz 
is nothing but the simplest PTA, the T® approximant. 

The purpose of this paper is twofold: first we want to ex- 
tract the slope and the curvature of the Transition Form Factor 
using a sequence of P^ approximants as fitting functions to 
the available experimental data. We demand an assignment of 
a systematic error to our predictions. Second we estimate the 
impact of our results on the light-by-light (LBL) contribution 
to the hadronic process on the muon g — 2. We also com- 
ment in passing on the recent proposal to measure the transi- 
tion form factor at low energy using the BESIII experiment. 

We proceed as follows: In section|II]we study the reliability 
of the PA method as fitting functions to proceed then on sec- 
tion|ni]to analyze the real data. In section|lV]we consider the 
impact of the previous result on the LbyL contribution to the 
muon g — 2. We finally collect all the results on the Conclu- 
sions section. 



II. TESTING THE METHOD WITH A MODEL 

Before applying the method to the experimental data to ex- 
tract the slope and the curvature of the TFF, we want to test its 
reliability with a particular model. Since it has not been pos- 
sible to describe rigourously the TFF from basic principals, 
several models have been developed during the last years with 
the purpose of analyzing the space-like data to extract fun- 
damental QCD properties. In these Refs. [16-38] we try to 
summarize the large effort done on this purpose. 



Considering this variety of models we examine instead of 
just one, three of them that we think are representative of the 
large amount of work done on this respect. Since our intent 
is to show the properties of our method, the selected models 
should describe well the experimental data but keep the com- 
plexity at a manageable level. This exercise will also provide 
a way to estimate the systematic error of our approximations. 

For easy of reading we comment here about the first model 
and relegate the other two to the Appendix [A] 

The first model considered is motivated by a Quark Model 
(e.g., ||T9H2~ni36l, see also fl28) for other log(g 2 /M 2 ) related 
models), although it can also be inspired by the lowest order 
perturbative QCD(pQCD) with a flat pion distribution ampli- 
tude (see, for example, GUI ) or even by the BABAR fitting 
function 0. We named the model "log-model": 

M 2 ( O 2 \ 

We 2 ) = ^^iog(i+|,j, o) 

with M 2 = 0.6GeV 2 and f % = 92MeV. 

Expanding F^yt^Q 2 ) in Eq. |3j in powers of Q 2 we obtain 

F*Pf r (Q 2 ) =a - ai Q 2 +a 2 Q 4 -a 3 Q 6 + G(Q % ) , (4) 

with known values for those a, coefficients (in particular ao = 
), as shown in the last column of Table |lj 

In order to illustrate the utility of the PA as a fitting func- 
tions we simulate the situation of the experimental data [iJH 
with the model by considering the function Eq. <[3j evaluated 
at 22 points in the region 0.7 < Q 2 < 5.5 GeV 2 , 16 points in 
the region 5.5 < Q 2 < 12.5 GeV 2 and 14 more points in the re- 
gion 12.5 < Q 2 < 35 GeV 2 . On top of these set of data points 
we add the value of 7^0^(0,0) = ■ All these data points 
have zero error because we want to obtain a pure systematic 
error on our fitting functions. 

We construct a sequence of Pi(Q 2 ) approximants with un- 
known coefficients as defined in Eq. (|2]) and then we fit the set 
of data which yields a predictions for the a, coefficients. The 
results are shown in Table [i] where we go up to the Pj 5 . The 
first Pf has only two parameters (ao and a\) and then a 2 is not 
a fitted but predicted through expansion. We also include on 
this table the position of the pole of each PA and the reader 
should notice how these poles, although showing a conver- 
gence pattern, differ from the lowest-lying vector mass used 
in the model of Eq. 

As expected fl2l . the sequence of PA converge to the exact 
result in a hierarchical way (much faster for ao than for a\ and 
so on), achieving with the last PA Pf a relative error of 0.04%, 
5.6% and 21.0% for ao,a\ and a 2 respectively. 

Similar results can be found by using as a fitting functions 
a sequence of PTAs as we said in the introduction. Thus, fix- 
ing the pole of the T(- at s p = M 2 = (0.77) 2 GeV 2 , we obtain 



1 The model of Eq. o3 1 has a branch cut starting at Q 2 < —M 2 . 
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P° 

1 


1 


P\ 

1 


1 


P} 

1 


1 


Fjfiyty (exact) 


a {GeV- 1 ) 


0.2556 


0.2694 


0.2734 


0.2746 


0.2751 


0.2752 


0.2753 


a x (GeV- y ) 


0.1290 


0.1716 


0.1935 


0.2051 


0.2124 


0.2166 


0.2294 


a 2 (GeV- 5 ) 


0.0651 


0.1147 


0.1492 


0.1725 


0.1898 


0.2013 


0.2549 




1.41 


1.22 


1.14 


1.09 


1.05 


1.03 


0.77 



TABLE I: aQ,a\ and a% low-energy coefficients of the log-model in Eq. |3j fitted with a P[{Q 2 ) and its exact values (last column). We also 
include the prediction for the pole of each P X (Q 2 ) (s p ) to be compared with the lowest-lying meson in the model. 



for the a relative error of 3%, 34% and 92% for ao,ai and 
ci2 respectively. These results could be easily improved if in- 
stead of fixing the pole of our PTA on the starting point of the 
branch cut, i.e., at s p — (0.77) 2 GeV 2 , we fix it at a different 
s p > (0.77) 2 GeV 2 . For example, if s p = lGeV 2 (value moti- 
vated by the result obtained with the previous P^) the relative 
errors turn out to be 0.15%, 2.3% and 14.4% for ao,a\ and 
ci2 respectively. Since the PTA's prediction are very similar to 
PA's ones, we do not show explicitly the corresponding table. 
This simple exercise shows that fixing the pole of our approx- 
imant to the physical resonance, as in the VMD case, might 
not be the best strategy to follow for low-energy constant pre- 
dictions, as extensively studied in Ref. l39l . 

The nice convergence pattern shown by our PA sequence 
should not be a surprise since it turns out that our model 
Eq. Q is a Stieltjes function and thus the convergence of the 
PA sequences is guaranteed by Pade Theory [40|. 

On the other hand, it has been recently considered in 
Ref. PD the possibility of KLOE-2 experiment at Frascati 
to measure the TFF at very low energies in the space-like re- 
gion (for 0.01 < Q 2 < O.lGeV 2 ) and the width r,j0^ yy at the 
per cent level. This new low-energy data may reduce our sys- 
tematic error for the PA P^ to 4.2% and 18% for a\ and a 2 re- 
spectively. An even better result might be obtained when the 
BES-III experiment at the e + e~ collider BEPC-II in Beijing 
will cover from the low-energy range up to CELLO energies, 
i.e., up to Q 2 ~ 0.7GeV 2 (which will turn out on systematic 
errors less than 3% and 15% for a\ and a 2 resp. considering 
the feasibility study for BES-III performed in 11421 ). Indeed, 
the 77 physics program at BES-III for the measurement of 
pseudoscalar TFFs will allow to cover a wide Q 2 range below 
lOGeV 2 , the gap between KLOE-2 and CLEO experiments. 

We analyze two more models in the Appendix [A] using the 
same technic explained here and the similar results obtained 
with the three of them give us a confidence on our fit proce- 
dure. This exercise allows us to assign a systematic error for 
each element on the PA and the PTA sequences. To ascribe a 
particular (and conservative) systematic error and taking into 
account we do not know the structure of the whole TFF, we 
select the worse of the three cases as a guidance. For PA Pj , 
5.6%, and 21% and for PTA T?, 5.4%, and 20%, as a relative 
systematic errors for a\ and a 2 respectively. 



III. FITS TO REAL DATA 

With all the tools developed so far we can now proceed to 
analyze the real TFF. For this purpose we use all the avail- 
able experimental data in the space-like region, which may be 
found in Refs. |Q]-0), and also the recent measurement of the 
r n o^yy decay width by the PrimEx Collaboration [43 1. 

The form factor for real photons is related to the n° — > yy 
decay width: 

F^yy{q\ = 0,q 2 = 0) = -±_T K 0^ r (5) 

with a = a em = 1/137.0356. 

The experimental world average collected in the PDG ta- 
bles ED is r™ G = 7.74 ± 0.48 eV, although we use here 
the PrimEx Collaboration result |43 1, which using a Primakoff 
effect experiment at JLab, has improved significantly the ac- 
curacy, reporting the value F K o^ n = 7.82 ± 0. 14 ± 0. 17 eV. 



A. Fits with the rational approximants 

The fits with the P^ sequence to the space-like data points 
in Refs. |Q}|4) determine those coefficients that best inter- 
polate them. As always, when fitting experimental data one 
should find a compromise between the increase of fit errors 
and decrease of systematic ones when increasing the order 
L of the f[ L . Figure [l] shows the experimental data obtained 
by CELLO (grey squares) 2 , CLEO (red diamonds), BABAR 
(blue circles) and Belle (brown triangles) together with the 
pQCD prediction (horizontal black dashed line). The red 
curve on FigjT|is our best approximant, the Pj\ 

In Fig. [2] we show the results for the prediction of the slope 
and curvature parameters a K and b K with the P^ up to L = 5. 
Approximants with L > 5 have the new coefficients compati- 
ble with zero and then do not introduce new information with 
respect to Pj 5 . The internal errors shown in Fig|5]are only sta- 
tistical and the external ones are a quadratic combination of 



2 CELLO data points Z), are extracted from Ref. 1 1 1 using the following nor- 
malization: Dj = ( aftffis ) ' 2 > with the N i = Fl( Gil'"" P rovided in that 
reference and a = 1/137.036. 
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FIG. 1: ThePf sequence compared with the y* y — > ;r° Transition Form Factor data from CELLO (grey squares) 1 1 1 , CLEO (red diamonds) 1 2 1 , 
BABAR (blue circles) |3| and Belle (brown triangles) |4): 



Zj (orange dashed), Pj (green dotted), Pf (brown short-dashed), (blue long- 



dashed), Pj (black dot-dashed) and Pj (red solid). Black dashed line indicates the pQCD result. 



statistical and systematic errors, the latter determined in the 
previous section. For completeness we also ascribe a 45% of 
systematic error to the PDG slope value 3 . The curvature pa- 
rameters has never been measured so for easy of comparison 
we expand the VMD fit used by the CELLO collaboration up 
to that order with the corresponding systematic error. 

As expected from the models studied, we see in these fig- 
ures a nice convergence pattern for both a K and b K . 

The PA Pf yields 

fl7C = 0.0340(35) stat (19) sys , (6) 

and 

^ = 1.20(28) stat (25) sys xKT 3 , (7) 

with a x 2 /d.o.f. — 0.80, where the systematic error is esti- 
mated from the previous section (5.6% for a n and 21% for 
bjt). We also extract the position of the PA pole s p = cil/cil+x. 
This ratio is shown in Fig (3] together with a band correspond- 
ing to the physical value M p ± r p /2 where M p = 0.7755GeV 
and F p = 0.155GeV is believed to be the dominant reso- 
nance contribution. For the Pf the pole is located at ^Jsp — 

0.75+^GeV, well in this band. 

It is interesting to notice the slightly larger results for the 
slope obtained with Py with L > 1, manifesting the need of 
a systematic procedure for going beyond VMD. It turns out 
that the larger L, larger the sensibility of the Pj L to the high- 
energy data. In this respect, the recent Belle data is crucial 
to obtain an accurate low-energy prediction since up to now 



3 Again, this systematic error is obtained comparing the VMD result with 
the exact one in TableslIIIVI andfv] 



Data a K (10 3 )fc g ^/g X 2 /d.o.f. 

All 0.0340(35) 1.20(28) 0.75+jJ;[^ 0.80 

CELLO+CLEO+BABAR 0.0348(39) 1.26(32) 0.73+oog 0.61 
CELLO+CLEO+Belle 0.0326(39) 1.08(30) 0.76+[}^ 0.49 

TABLE II: Slope and curvature of the TFF predictions with a Pj with 
different sets of data. 



BABAR data was dominating the high-energy region (see, for 
example, Ref. lfL3l for a preliminary study without Belle data). 

Our final result is a fit to all the available data but for a 
deeper understanding of PA as a fitting functions for high- 
energy data we consider two different scenarios where first no 
Belle data is considered and second no BABAR data is con- 
sidered (Table [III). Surprisingly enough, the results shown in 
Table [ill (where -A" stands for all the available data) are nicely 
compatible within errors although with slightly different cen- 
tral values. All the results are obtained with a Pj\ 

Finally, we want to apply a last test of robustness to the 
method which is fits to subsets of data should return compati- 
ble results, i.e., fitting data up to lOGeV 2 , up to 20GeV 2 and 
up to 36GeV 2 (all the data) should be the same (unless there 
are unknown problems with the data like normalization or sys- 
tematics). The results are shown in Table [HT] where we also in- 
dicate the best P\" that fits the particular subset of data. These 
results are nicely compatibles, otherwise we should take the 
difference as a new source of systematic error. 

For illustrative purpose we show on Fig.|4]the result when 
fitting data up to lOGeV 2 . In this case we include also the fea- 
sibility study for BES-III experiment performed in Ref. |42l . 
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FIG. 2: a K (left) and i^(right) predictions with the P\ up to L = 5. The internal band is the statistical error from the fit and the external one is 
the combination of statistical and systematic errors determined in the previous section. 
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FIG. 3: Position of the pole ^fsj, for the different Pf. For comparison 
we also show (gray band) the range M p ± Tp jl corresponding to the 
physical p— meson value. 
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FIG. 4: The Pf sequence compared with the Y*y^r 7t° Transition 
Form Factor data up to lOGeV 2 GHH 02): Pf (orange dashed), 
P\ (green dotted), P 2 (brown short-dashed), Pj 5 (blue solid). Black 
dashed line indicates the pQCD result. 





best PA 




X 2 /d.o.f. 


Data up to lOGeV 2 


f 


0.0364(51) 


0.53 


Data up to 20GeV 2 


Ft 


0.0327(35) 


0.69 


All 


p l 


0.0340(35) 


0.80 



TABLE III: Slope of the TFF prediction with with different sets of 
data as described in the main text. 



With a we obtain 4 a n = 0.036(6) and b n = 1 .41 (65) where 
the errors are statistical and systematical with % 2 jd.o.f . = 
0.53, to be compared with the results on Table [H] In this sce- 
nario the pole of the P\ is located at ^/sj, = Q.Ty^ '^GeV. 



PA with larger L do not introduce new information. BES-III data will be 
crucial to improve on this result. 



B. Other Pade Approximants 

1. Pk Pade Approximants 

The experimental data so far considered range up to 
36GeV 2 then a natural extension of the previous analysis 
would include higher resonances although the form factor is 
believed to be dominated by the p (770) meson. In such a way, 
the consideration of two-pole p£ will give us a way to asses 
any possible systematic bias in our P^ analysis. 

In this case our best appro ximant is the This approxi- 
mant yields 

a* = 0.0324(20) and b„ = 1.07(15) X 10~ 3 , (8) 

with a x 2 jd.o.f . = 0.71, nicely compatible with our previous 
determination in Eqs. ([6]) and Q. Despite this result, the poles 
of that approximant are located at s p \ = 0.53(6) — /0.01(1) 
and s p 2 = 0.56(2) + /0.01(1), where we can see a certain pa- 
rameter space region where the poles may became eventually 
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complex-conjugated 5 . 



C. Final Result 



2. 7j Pade-Type Approximants 



The results shown in Eqs. (|6fTTj) agree quite well and com- 
bining them our final weighted average result yields: 



On the other hand, since the value of the physical p —meson 
mass is well known, it is natural to attempt to include this 
information in our analysis through the PTAs. We have seen 
on the previous section, however, that locating the pole of a 
PTA exactly at the physical counterpart is not the best strategy. 
In fact, we learnt that s p > M 2 but it is not clear what particular 
value we should use. To evaluate a possible systematic error 
on this choice, we range the PTA pole in between the band 
drawn by the PA results in Figj3j i.e, ^/Sp = 0.73 - 0.83GeV. 

With a T\ sequence we go up to and obtain also a nice 
and smooth convergence pattern for both a K and b n parame- 
ters. With our best PTA, we obtain 



a K = 0.0302(28) and b K = 0.92(18) x 10" 



(9) 



with a x 2 /d.o.f. = 0.78 — 0.87, where the errors are mainly 
the systematics of the pole range. 



'N+i 



Pade Approximants 



^ = 0.0324(12)^(19) 



sys > 



and 



^ = 1.06(9W(25),„xlO- 



(12) 



(13) 



to be compared with other theoretical determinations: from 
a Regge analysis, a n = 0.032(1) HTJ; from ChPT at the 
loop level with pi = M p , a K = 0.036 PHTl : from a study of 
the Dalitz decay 7t° -> e+e~y, a n = 0.029(5) B31 : from a 



hard- wall holographic models of QCD, a K « 0.031 114611 and 
a K ss 0.035 [47 1; from a soft-wall holographic model of QCD, 
a n = 0.024(5) J48] 6 ; and finally from the compilation of holo- 
graphic models on Ref. [49], a K = 0.031(6), where the error is 
estimated by the spread of the different results obtained from 
these models. 

In the next section we explore possible consequences of our 
final results in Eqs. ( 12 1 and ( 13 1 on the light-by-light scatter- 



ing contribution to the anomalous magnetic moment of the 
muon. 



As suggested in the Introduction, we may attempt to include 
the asymptotic behavior of the TFF 1 14 1 in our fits by consid- 
ering a P$ + i sequence. With these approximants we can go 
up to the P 2 which yields 

a n =0.0331(45) and b % = 1.11(27) x 10~ 3 , (10) 

with a x 2 /d.o.f. =0.73. 

The P 2 approximant has the right fall-off as Q~~ 2 but the 
corresponding coefficient (which reads 0.17 ± 1.8GeV) is not 
correctly predicted due to its large statistical error. This 
asymptotic coefficient is known to be 2f K by first princi- 
ples fl4l . It seems logical to trying to include this information 
on the P 2 . Using the asymptotic coefficient when construct- 
ing the P 2 , we obtain a constrained approximant called P' 2 
which after fitting the TFF data yields, for the low-energy co- 
efficients, 



= 0.0332(25) and b x = 1.13(19) X 10" 



(11) 



with a x 2 /d.o.f. =0.70. 

It is remarkable that these results with the P' 2 , which makes 
use of all the experimental data and the asymptotic limit at 
once, are nicely compatible with all the previous results. This 
approximant seems to suggest that the scale where the pQCD 
should be applied is much further away than the last BABAR 
and Belle data points. 



IV. IMPLICATIONS ON THE HADRONIC 
LIGHT-BY-LIGHT CONTRIBUTION TO THE (g 2) M 



We can use the results in Eqs. ( 12 1 and ( 13 1 to constrain any 



model that estimates the pion-exchange piece to the Light-by- 

Light scattering contribution to the (g — 2) ll , the ' ' term. 
As an example, we consider the so called LMD + V model 
(defined in Ref. 11501 ) to account for that contribution: 



F LMD+V ln 2 n 2\ _ 

h -Q]Q 2 2 {Q] + Q 2 2 ) + hi(Q] + Q\) 2 + h 2 Q 2 Ql - h 5 (Q 2 + Q 2 ) + in 



(Q 2 +M 2 )(Q 2 



-M 2 ;)(Q 2 +M 2 )(Q 2 - 



(14) 



The TFF is related to the LMD+V model Eq. ([T4} when 
one of the photons on the latter is on-shell. That means 
we cannot fix all the free parameters (hi, with i — 1,2,5,7, 
and M Vl ,M Vl ) on this LMD+V model at once. We need 
more information, for example, from the high-energy region 
(Q 2 F n o rr (Q 2 ,0) = 2f n , Ref. d). If we match the high- 
energy limit, we find hi = and h$ = —6My i My 2 . The axial 
anomaly on the low-energy limit fixes /17 = — - N £ 2 M v My^ . 



' Complex-conjugate to render the approximant real. 



6 This number is obtained through the large-A^ limit relation = ft °^^ 2 

l45l with C^ 2 = 6.3 x 10 -3 obtained in [48]. Indeed, with our final value 
for a K we predict = 8.4(9) x l(r 3 GeV- 2 . 



7 



With these results and hi — as suggested in Ref. 11501 . we 
can use the slope and the curvature of the TFF to fix My l and 
M Vl . WefindM 2 ,, = 0.33(1 l)GeV 2 and M Vi = 0.94^GeV 2 



'Vi 

'H 1 



and w e obtain a^ byL ' K = 5.4(5) x 10 



10 



saturated by infinitely many sharp meson states. In the par- 
ticular case below, the pion couples first to a pair of vector 
mesons V p and V m which then transform into photons. Thus 
we have: 



V. CONCLUSIONS 

In this paper, we analyzed the collection of all the 
experimental data on the n°y*Y Transition Form Factor 
at low-energies with a model-independent approach based 
on Pade Approximants and we obtain the slope a K = 
0.0324(12)^(19)^ and curvature b K = 1.06(9)^(25)^ x 
10~ 3 of the form factor. The method is simple and system- 
atic and provides a model-independent estimation of all the 
systematic errors. We analyzed the impact at low-energy of 
the Belle and BABAR high-energy data and also the future 
BES-III data. We also evaluate the implications of these re- 
sults on the pion-exchange contribution on the light-by-light 
scattering part of the anomalous magnetic moment of the 
muon. Using the well-known LMD+V parametrization and 
the Pade Theory technics we estimate that contribution to be 

o^=5.4(5)xl0- 10 



E 



v y a {l\-^M~ M vJ 



+ {q\^q%), (Al) 



where Fy and Fy m are the current- vector meson couplings and 
G n VpV a is the coupling of two vector meson to the pion. The 
dependence on the resonance excitation number n is the fol- 
lowing 

M Vp = Ml m = M 2 + nA 2 , and F Vp = N L V a = F . (A2) 

The combination of sums in Eq. ( |A1[ > can be expressed in 
terms of the Digamma function \j/(z) — 4-logT{z): 



F n o rr (q u q 2 ) = F n o rr (Q ,A) = 

c r (m 2 e 2 (i+A) \ 

N t AQ 2 y\A 2 2 A 2 J 



(A3) 



(M 2 



Q 2 (l-A) 
2A 2 
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Appendix A 



For completeness we also studied two more models for the 
F J fiy t y,{q\,q2). The first is based on the Regge theory and the 
second on the light-front holographic QCD. After generating 
a set of zero error data points for each model we fit the data 
withPj L (g 2 ) and T^Q 2 ) sequences. 



1. Regge-model 

We consider first a Regge model based on the laige-N c 
limit, N c been the number of colors, (see, for example, 
|[T6l [T7l l30l l35l were similar large-A^ c models are used to fit 
directly the available data). In this limit, the vacuum sector 
of QCD becomes a theory of infinitely many non-interacting 
mesons and the propagators of the hadronic amplitudes are 



where Q~ = — ( 



-qi, 



,A = 



2 2 



and c a constant. 



To reassemble the physical case we consider N c = 3, A 2 = 
1 .3GeV 2 (as suggested by the recent light non-strange qq me- 
son spectrum analysis OP ). A = 1 (which means q 2 = 0), 
M 2 = (0.8) 2 GeV 2 and the constant c in such a way that the 
anomaly 7^0^(0,0) = 4 J,^. is recovered. 

Eqs. \A\\ and ( |A3| > use the large-A^ : and chiral limits and 
thus have an analytic structure in the complex momentum 
plane which consists of an infinity of isolated poles but no 
brunch-cut (as does have the Log-model of section [EIJ, i.e. 
they become meromorphic functions. As such, they have a 
well-defined series expansion in powers of momentum around 
the origin with a finite radius of convergence given by the first 
resonance mass. It is well-known [52 1 and largely explored 
in the context of Large-A^ ||39l l53l than the convergence of 
any near diagonal PA sequence to the original function for 
any finite momentum, over the whole complex plane (except 
perhaps in a zero-area set) is guaranteed. 

For meromorphic functions such as Eqs. ( A1|A3 i, another 
important result of Pade Theory applies here, the Montessus 
de Ballore's theorem lfl3l [151 . which states that given a cer- 
tain analytic function f(z) at the origin which is meromorphic 
with exact M poles in a certain disk on the complex plain, the 
sequence of PA converges uniformly to f(z). In practice, pro- 
vided M is known (M = 1 in our case), the Montessus' theo- 
rem asserts convergence for the sequence of M-pole Pade Ap- 
proximants Pf^. These convergence theorems are confirmed 
by the good results collected on Table IV where after generat- 



ing a set of zero-error data points with the model of Eq. ( A3 1, 
we fit them with the PA sequence and obtain the predictions 
for the a, coefficients. 
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P? 


Pi 


P? 


n 


P\ 


P\ 


fjOyi, (exact) 


a a {GeV~ l ) 


0.2672 


0.2730 


0.2746 


0.2751 


0.2752 


0.2753 


0.2753 


ui(G«V- 3 ) 


0.2662 


0.3121 


0.3338 


0.3457 


0.3529 


0.3571 


0.3678 


a 2 (GeV- 5 ) 


0.2652 


0.3600 


0.4244 


0.4616 


0.4868 


0.5030 


0.5550 




1.00 


0.92 


0.87 


0.86 


0.85 


0.84 


0.80 



TABL E IV: a$, a\ and a 2 low-energy coefficients of the Regge-model 
in Eq. i A3 i fitted with a P\[Q 2 ) and its exact values (last column). 
We also include the prediction for the pole of each P\{Q 2 ) (s p ) to be 
compared with the lowest-lying meson in the model. 





P? 


Pi 


P\ 


Pi 


1 


P? 


fyy r (exact) 


a (GeV- 1 ) 


0.2791 


0.2774 


Q.21M 


0.2759 


0.2756 


0.2754 


0.2753 


ai(GeV~ 3 ) 


0.3571 


0.3362 


0.3213 


0.3108 


0.3033 


0.2986 


0.2856 


a 2 (GeV- 5 ) 


0.4567 


0.4031 


0.3643 


0.3358 


0.3148 


0.3009 


0.2535 


sfr(GeV) 


0.88 


0.91 


0.94 


0.96 


0.98 


1.00 


1.16 



TABLE V: an, a] an d at low-energy coefficients of the Holographic- 
model in Eq. iA4i fitted with a P^Q 2 ) and its exact values (last 



column). We also include the prediction for the pole of each Pf'iQ 2 ) 
(s p ) to be compared with the lowest-lying meson in the model. 



With the PA P{ we obtain a relative error of 0.02%, 2.9%, 
and 9.4% for ao,fli and ci2 respectively. The inclusion of the 
feasibility study at BES-III ll42l decreases the error down to 
2.4% and 7.9% for a\ and a 2 resp. With a PTA sequence 
the results return 0.02%, 0.7%, and 0.8% for ao,ai and a 2 
respectively when the PTA pole is located at s p = 0.70GeV 2 



2. Holographic-model 

Finally, as a third model be analyze a simple holographic 
confining model presented in ||32l (and also explored in 
Refs. 0JD E2 HIES), based on light-front holographic QCD 
where the correct small Q 2 behavior (in order to simulate con- 

7 We do not consider higher-twist components to keep the model easy to use. 



finement) is introduced using the dressed current (see 11321 for 
details) 7 . 

In this context, the TFF is defined as 



(Q 2 ) 



X 2 fnJo (1 



dx 



(A4) 



Fjfiyy 



where P q q is the probability of finding the qq component in 
the pion light-front wave function. To reproduce the anomaly 
(0) = \/{An 2 f K ), we impose P qi = 0.5. 
This model reproduces quite well the transition form factor 
data up to lOGeV 2 but disagrees in particular with BABAR's 
large Q 2 data (although compatible with Belle data within 
errors), specially because the model is reaching its asymp- 
totic prediction (Q 2 F K o rr (Q 2 -> 00) = 2f K OH) already at this 

medium- Q 2 region. Another interesting feature of this model 
is that no convergence theorem from Pade Theory is known 
for this kind of function and then the Pade convergence is 
not guaranteed in advance (in contrast to the previous Regge- 
model). It represents a robustness test of our method. 

After generating again a set of zero-error data points with 



the model of Eq. ( A4 1, we use the PA sequence to fit these 
data and to obtain again the predictions for the a, coefficients. 
We collect the results in Table IVl 



With the PA Pf 



we obtain 0.04%, 4.6%, and 18.7% as 
a relative errors for oq, a\ and a 2 respectively. With the 
inclusion of the feasibility study at BES-III |42| . we go 
down to 4.3%, and 17.1% for a u a 2 resp. With the PTA 
sequence (the approximant pole located at s p = lGeV 2 ) we 
obtain 0.6%, 4.8%, and 19.2%, respectively. Although no 



convergence theorem for this kind of function in Eq. (A4i is 
known, the convergence of our PA sequence is clear. That is 
one of the most interesting features of the PA methods which 
is the convergence may occurs beyond expectations. 
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